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LISA Pathfinder is a science and technology demonstrator of the European Space Agency within 
the framework of its LISA mission, the latter aiming to be the first space-borne gravitational wave 
observatory. The payload of LISA Pathfinder is the so-called LISA Technology Package, which 
is designed to measure relative accelerations between two test masses in nominal free fall. The 
diagnostics subsystem consists of several modules, one of which is the magnetic diagnostics unit. 
Its main function is the assessment of the differential acceleration noise between the test masses 
due to magnetic effects. This subsystem is composed of two onboard coils intended to produce 
controlled magnetic fields at the location of the test masses. These magnetic fields couple with the 
remanent magnetic moment and susceptibility and produce forces and torques on the test masses. 
These, in turn, produce kinematic excursions of the test masses which are sensed by the onboard 
interferometer. We prove that adequately processing these excursions, the magnetic properties of 
the test masses can be estimated using classical multi-parameter estimation techniques. Moreover, 
we show that special processing procedures to minimize the effect of the multi channel cross-talks 
are needed. Finally, we demonstrate that the quality of our estimates is frequency dependent. We 
also suggest that using a multiple frequency experiment the global estimate can be obtained in such 
a way that the results of the magnetic experiment are more reliable. Finally, using our procedure 
we compute the the contribution of the magnetic noise to the total proof-mass acceleration noise. 

PACS numbers: 04.80.Nn, 04.30.-w, 07.87.-|-v, 06.30.Ka, 07.05.Fb 



I. INTRODUCTION 

LISA Pathfinder (LPF) is a science and technol- 
ogy demonstrator programmed by the European Space 
Agency (ESA) within its LISA mission activities [T]. 
LISA (Laser Interferometer Space Antenna) is a joint 
ESA-NASA mission which will be the first low frequency 
(milli-Hz) gravitational wave detector, and also the first 
space-borne gravitational wave observatory. The payload 
of LPF, the LISA Technology Package (LTP), will be 
the highest sensitivity geodesic explorer flown to date. 
The LTP is designed to measure relative accelerations 
between two test masses in nominal free fall (geodesic 
motion) with a differential acceleration noise budget 
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in the frequency band between 1 niHz and 30 niHz |21 . 

Magnetic noise in the LTP is allowed to be a signifi- 
cant fraction of the total mission acceleration noise: up to 
1.2 X 10"^'* ms^^ Hz^i/^ is apportioned to magnetic ef- 



fects, i.e., 40% of the total noise, 3x10 
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see Eq. ^. This noise occurs because the residual mag- 
netization and susceptibility of the test masses couple to 
the surrounding magnetic fleld, giving rise to a fluctuat- 
ing force which is given by: 
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in each of the test masses. In this expression B is the 
magnetic field in the test mass, x ^-nd M are its magnetic 
susceptibility and residual density of magnetic moment, 
respectively, and V is the volume of the test mass, /ig is 
the vacuum magnetic constant, 47r x 10~^ mkgs~^ A~^, 
and (■ • • ) indicates test mass volume average of the en- 
closed quantity. Finally, (5B represents the fluctuation 
of the magnetic field, and (5V stands for the fluctuation 
of the gradient [1]. Quantitative assessment of the mag- 
netic noise in the LTP, i.e. evaluation of Eq. ([2| clearly 
requires a real-time monitoring of the magnetic field and 
an accurate knowledge of the magnetic characteristics of 
the test masses. 

The determination of the magnetic characteristics of 
the test masses (remanent magnetic moment and sus- 
ceptibility) must be done in flight because their mag- 
netic properties may change due to launch stresses and 
other circumstances. This will be done injecting con- 
trolled sinusoidal magnetic fields at the positions of the 
test masses and appropriately processing the resulting 
kinematics, which will be obtained from the readings de- 
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livered by the onboard interferometer. Although the ba- 
sic design of the magnetic experiment is well settled [5] , 
due to the high complexity of the LTP experiment, more 
in-depth analyses based on a more realistic modeling are 
necessary to assess its feasibility and performance. The 
purpose of this paper is, precisely, to fill this gap. In 
particular, we model in a realistic way the kinematics of 
the test masses and we evaluate the expected quality of 
the estimates of the magnetic moment and susceptibility. 
Specifically, we take into account several effects — like 
the cross-talks between some of the channels of the in- 
strument, or the frequency-dependent parameters of the 
control loops governing the dynamics of the test masses 
— that previous analyses disregarded. All these effects 
depend on the frequency used to excite the test masses. 
Hence, the quality of the estimates of the magnetic data 
depends sensitively on the excitation frequency, since 
the satellite does not behave equally across the complete 
measurement bandwidth. Therefore, it is important to 
determine the quality of the estimates across the com- 
plete measurement bandwidth, and the frequencies that 
deliver the best estimate of the magnetic parameters. 

The paper is organized as follows. In Sect. |Tl]we give 
a brief description of the magnetic experiment intended 
to estimate the magnetic properties of the test masses. 
Then, in Sect. |IIl] we briefly present the dynamical model 
of the satellite. It follows Sect. IV where we discuss the 
estimation model and the estimation procedures used in 
this work. In Sect. [V] we present the sensitivity of this 
model to different hardware systems of the satellite and 



in Sect. |Vl] we evaluate the frequency dependence of the 
experiment, and we optimize its performance with re- 
spect to the excitation frequency. Sect. |VlI| is devoted to 
analyze the robustness of our results. Sect. |VlTl] deter- 
mines the expected accuracy of the magnetic contribu- 
tion to the total proof-mass acceleration noise. Finally, 
in Sect. [IX| we summarize our main findings, we discuss 
the significance of our results, and wc draw our conclu- 
sions. 



II. EXPERIMENT OVERVIEW 

The two test masses are located at the center of each 
inertial sensor — the two towers in Fig. [T] — and are the 
end mirrors of the Optical Metrology System, that senses 
the positions and attitudes of the test masses. The op- 
tical bench of the interferometer can be seen in Fig. [T] 
as well. In fact, one of the test masses is the reference 
free floating body to perform the translation and atti- 
tude control of the spacecraft. The x-axis of the experi- 
ment is the axis connecting the two test masses centers, 
and it goes from test mass 1 to test mass 2. The 2-axis 
points towards the solar panel (parallel to the two iner- 
tial sensor towers and upwards in Fig.[T]) and, finally, the 
y-axis closes the right-handed Cartesian coordinate sys- 
tem. The test masses are made of an alloy of Pt (27%) 
and Au (73%), their dimensions are 46 x 46 x 46 mm and 
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FIG. 1: A schematic view of the payload of LISA Pathfinder, 
the LTP. The Inertial Sensors (two vertical towers) host the 
two test masses. The four floating boxes correspond to the tri- 
axial fluxgate magnetometers, and the two induction coils are 
placed next to each of the test masses. The optical bench of 
the interferometer is located on the horizontal plane between 
the Inertial Sensors. 



their weight is 1.95 kg. To comply with the top science 
requirements, the test masses must have certain proper- 
ties. For the purpose of the present work the two most 
important properties are the remanent magnetic moment 
and the susceptibility. The remanent magnetic moment 
must be |m| < 2.0 x 10"^ A m^. Since the volume of 
the test masses is V — 0.046'^ m"^, the density of mag- 
netic moment must be then |M| < 9.451 • 10"'' A/m. 
The susceptibility of the test mass can be suitably repre- 
sented by a complex number, X = Xo + ^Xc: where Xo is 
its real component and Xo(w) is a frequency-dependent 
imaginary term which is due to the eddy currents on the 
test mass I3| . The requirement on the value of the real 
component is Xo < 2.5 x 10"^. 

As mentioned, to measure the remanent magnetic mo- 
ment and the susceptibility of the test masses a controlled 
magnetic field will be injected at the position of the test 
masses. This magnetic field produces forces and torques 
which excite the kinematics of the test masses. Studying 
the motion of the test masses, namely their displacement 
and rotation, allows to estimate the three components of 
the magnetic moment and the susceptibilities of the test 
masses. 



A. The injected magnetic fields 

The magnetic field at the position of the test masses 
will be generated by the injection of sinusoidal currents to 
the onboard coils. These onboard coils are placed next 
to each of the inertial sensors towers, see again Fig. [l] 
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The two circular induction coils are made of a titanium 
alloy (TigAUV), and have N = 2 400 windings of radius 
r = 56.5 mm [3]. They are placed 85.5 mm away from 
the center of the respective test mass. The onboard coils 
are aligned with the x-axis of the test masses, thus, the 
magnetic field within the volume of the test masses has 
axial symmetry. Given a current fed to the coils I{t) = 
/q sin Wot, the resulting magnetic field (and its gradient) 
will oscillate at the same frequency. Therefore we write. 
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The field produced by the coils at the center of the test 
mass is 4.47 f/T, whereas the maximum environmental 
magnetic field expected during science operation is less 
than 100 nT. On the other hand, the magnetic field gradi- 
ent along the x-axis produced by the coils is 109.2 /iT/m, 
while the maximum magnetic field gradient required by 
the mission science specification is —5 f/T/m. Therefore, 
it is a safe assumption to neglect the effects of the envi- 
ronmental magnetic field with respect to the applied field 
by the coils. Thus, the forces and torques exerted on the 
test masses are computed as (B]: 
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where Bapp is the field produced by the coils, and r is the 
position vector that has the test mass center as origin. 

The term Re | ^"^p'^" Bo ie~*""*| is the induced magne- 
tization by the externally applied field. Note that the 
forces and torques depend on M, Xo and Xc- 

Considering Eqs. ([s]), Q and the x-component of 
the force acting on the test mass is 
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where we have used that sin^ wot = (1 — cos2wot)/2 and 
the 7r/2 rad phase due to the complex component of the 
susceptibility has been added as an argument in the cor- 
responding cos term. As can be seen from this equation, 
the linear acceleration of the test masses along the x-axis 



has two separate frequencies, one at ojq and the other 
at 2wo, and also a DC component. The 2wo component 
presents an in-phase component proportional to Xo and a 
quadrature component proportional to Xc- Particularly, 
the cjo component can be more explicitly written as: 
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where M^, My, are the components of the density 
of the remanent magnetic moment. If the test mass is 
homogeneous we have the simplified expression 
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since the y and z components of Vi?o,2: average to zero 
due to symmetry of the field of the coil. This leads to a 
force component along the x-axis that only depends on 
Mx, Xo and Xc- 

On the other hand, the torque acting on the test mass 
also has a similar behavior: 

N = (M X Bo + r X [(M-V)Bo]) y sinwot (10) 

+ ( r X ^ [Bo • V]Bo]\ V sin^ coot 
A*o / 
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In this case, it must be noted that, because of the symme- 
try of the applied magnetic field, the terms multiplying 
sin^ Wot and sinwot coswot in Eq. (10) vanish. The two 



rotation excursions detected by the interferometer using 
wavefront sensing are the rotations about the y-axis and 
z-axis. The magnitude of the rotation about the x-axis 
is smaller, and cannot be detected by the interferome- 
ter because the axis of rotation is aligned with the laser 
beam. Taking this into account, the two relevant torques 
for the experiment are: 
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These equations can be further simplified in the case 
of a homogeneous test mass. In this case, due to the 
axial symmetry of the magnetic field, the terms (i?o,z) in 
Eq. (11) and {Bo^y) in Eq. (12) vanish. Moreover, for the 
same reason the terms 
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in Eq. ( 11 ) also vanish, as do the terms 
dBo^y 



dx 



dBo^x 
dx 



dz 
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in Eq. (12). In these terms x, y and z are the three 



components of r. Hence, the torque about the y-axis 
only depends on and the torque about the z-axis 
only depends on My: 



N„ = ( Bo,, + z - X ^1^ )1^ sincjot (13) 
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Finally, we can cast Eqs. ([t]), ([T3| and (14) in the form: 



Fx = Xo fxuc + Mx fxi^g + Xo fx2^ + Xc fx 

Ny 
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where /ajj^^ is a constant function, fxi^g , and Hz^^ 

oscillate at loq and fx2^ and oscillate at 2aJo- 



III. DYNAMIC MODEL 



2. The sensing mechanisms (S) onboard LPF are the 
star tracker, the inertial sensors, and the interfer- 
ometer. Of special interest is the interferometer, 
which measures the distance between test mass 
1 and the spacecraft and between the two test 
masses [TO] . 

3. The controller blocks (C) are in charge of calculat- 
ing the appropriate commands to correct the po- 
sitions of the test masses and the attitude of the 
spacecraft. In science mode, there are two main 
control loops applied by the instrument. The first 
one — the drag free loop — takes the absolute mea- 
surement of the distance between test mass 1 and 
the spacecraft as a reference. It then calculates 
which forces should be applied to the spacecraft in 
order to counteract all disturbances and recreate a 
drag free environment for test mass 1. The second 
loop — the low frequency loop — takes as a ref- 
erence the differential measurement between both 
test masses and acts on the second test mass to 
avoid its collision with the spacecraft walls • The 
controllers have been designed to deliver very sensi- 
tive readings of the differential motion of both test 
masses between 1 mHz and 30mHz, the measure- 
ment bandwidth of the LPF mission f5] . These two 
control loops are implemented inside the onboard 
computer of the LPF. 

4. The actuators (A) are the physical systems that 
apply these commands to the test masses and to 
the spacecraft. The two actuator mechanisms ex- 
isting in LPF are the satellite micropropulsion sys- 
tem, which is composed by 12 micro-newton FEEP 
thrusters (Field Emission Electric Propulsion), and 
the capacitive actuators which consist of a set of 
electrodes that surround the test masses and exert 
controlled forces on them. 



The LTP instrument will react to the injection of the 
aforementioned forces and torques inflicted upon the test 
masses. This will result in specific kinematic excursions 
in both test masses. These kinematic excursions will de- 
pend on the instrument dynamics and will be sensed by 
the onboard interferometer. The LTP is a very complex 
instrument and its modeling has been presented in sev- 
eral references It can be modeled by splitting it in 
four main subsystems which are: 

1. The dynamical model (D) represents the evolution 
of the kinematic excursions of the two test masses 
placed inside the LTP and the kinematics of the 
spacecraft. This model takes into account the cou- 
pling of the motion of each of the test masses with 
the motion of the spacecraft and outputs the evolu- 
tion of the 15 degrees of freedom of the instrument 
(6 for each of the test masses and 3 more for the 
spacecraft) . 



This subsystem division is schematically shown in the 
block diagram of Fig. [2j For more detailed information 
of the system, the reader is referred to Refs. [THS]. 

In the magnetic experiment the input signals are the 
magnetic forces and torques (fmag), and the outputs are 
the readings of the interferometer (o). Therefore, using 
the block diagram scheme shown in Fig. ^ , we calculate 
the transfer function, which results in: 



o = H f , 



where 
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This transfer function depends on all the above described 
subsystems and represents the dynamical response of the 
instrument to the specific injected signals. 
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FIG. 2: Control system architecture of LISA Pathfinder. 
D stands for the dynamical matrix, S represents the sens- 
ing matrix of the interferometer, i.e. the matrix translating 
the position of the test mass, x, into the interferometer read- 
out, o (on stands for the readout noise). A represents the 
physics of the FEEP and the electrostatic actuators, and fi- 
nally C, is the controller matrix, implementing the drag free 
and low-frequency control loops. Oi represents the displace- 
ment guidance signals, an are the actuators noise and fa are 
the output forces and torques of the actuators, fmag are the 
magnetic forces and torques induced by the coils and fn are 
the environment force and torque noises disturbing the space- 
craft. 



IV. ESTIMATION MODEL 
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where the 3x3 matrix H is the transformation matrix 
from force/torque to displacement/rotation that repre- 
sents the closed loop dynamics of the instrument — see 
Eq. (17). This matrix is not diagonal, as it is assumed 



in the model in which the cross-talks are neglected — 
namely, Eq. (18 1. For instance, the effect of the torque 
about the y-axis and the z-axis on the Oxi2 displacement 
channel is relevant, and thus non-zero transfer functions 
HNy~¥xi2 and Hpf^^xi2 need to be considered. Hence, to 
estimate M^, My, M^, Xo and Xo, these transfer functions 
have to be known. This model is still a simplification, be- 
cause we do not include all the degrees of freedom, but 
it is certainly more realistic than that of Eq. ( 18 ), which 
is strictly one-dimensional. 



The estimation of the magnetic characteristics is per- 
formed processing the interferometer readings. To do so, 
we use the displacement of the differential channel {ox-^^^)^ 
the rotation about the y-axis (0,,^ ) and the rotation about 
the z-axis (o^J. If cross talks are disregarded, the read- 
ing of the displacement channel stems only by the effect 
of the magnetic force acting along the a;-axis. Fx- Anal- 
ogously, something similar occurs for the two torques in 
each of their respective axis. Thus, we can write: 



Oxi2 = Mx dx,^^ + Xo dx2^^ + Xc d^^^^ 
Om = ry,^^ (18) 
0^1 



where dx-^^^ , dx2^g and d^^^ are the respective transfor- 
mations from force to displacement of the signals fxi^„ , 

fx2^r, ^^'^ fx 



in Eq. (15), and analogously for 
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and rir, . Nevertheless, 



and r 

because of the high complexity of the LTP instrument, 
this model is not sufficiently realistic. In particular, it 
turns out that the cross-talks cause important biases in 
the parameter estimates. This is because the effect of 
the x-force in the rotation readings and the effect of the 
torques in the x-axis readings are not negligible. As a 
consequence, we used the full three-dimensional model 
of the experiment: 



A. Estimation procedure and bias correction 

The estimation procedure has been already described 
in Ref. [S]. However, in this paper we present an impor- 
tant modification to correct for the biases introduced by 
cross-talks. The full three-dimensional estimation model 



given by Eq. ( 19 ) may be regrouped as: 
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(20) 



where a are the cross-talks of the system (the matrix 
elements of H evaluated at the excitation frequency) . If 
we introduce primed quantities M^, My, and M^, then 
Eq. ( 20 ) can be written as: 



0x12 = M^ dx,^^ + Xo dx2^g + Xc ^2^0 

Or, = m', ry,^^ 
= My r^,^^ 



(21) 



and we estimate the values of AI^, My, M^, Xo and x^ 
applying standard single output least square techniques 



[5l[T6]. These values of M^, My and M^ are biased, and 
do not correspond to the true magnetic moment compo- 
nents, Mx, My and M^. Nevertheless, these biases can 
be corrected because we know the relation between them: 
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(22) 



Note that Eq. ( 19 1 provides the vahies of the elements 
of this matrix, and that the matrix is invertible. Ad- 
ditionahy, it is worth emphasizing that we only correct 
the components of the magnetic moment and no correc- 
tion is considered for the susceptibility (xo and Xc)- This 
is because the magnetic susceptibility is not affected by 
any cross-talk. It turns out that the previously outlined 
procedure corrects biases of around 1% in each of the 
magnetic parameters, which are sizable. Finally, we also 
mention that during the lifetime of the mission some of 
the telemetry channels may fail. Thus, it is important to 
know beforehand that single channel estimation is still 
possible and that it introduces biases of ^ 1%. On the 
other hand, by direct usage of the LTPDA toolbox, and 
in order to avoid additional estimation biases induced by 
the low frequency behavior of the instrument, we whiten 
the data and we eliminate its transients. These tech- 
niques are expected to be used in other experiments of 
the mission and are described in detail elsewhere [H [12]- 
[W]. 



wi = wi„o^,(l±6^J 

W2 = W2„o^,(l±6^J 

Sl2 = (5i2nom(1±^<5i2) 

GfeEP = GfeEPnom(1 =•= 

GcA = G'cAnom(1 ± ^Gca) 

Clearly, the effects of these uncertainties on the esti- 
mation of the magnetic parameters need to be assessed. 
To this end, for each of the nine transfer functions of H, 
we have computed the effect of the uncertainties on each 
of the parameters of the system. We have done this anal- 
ysis for values of b ranging from —0.2 to 0.2, and we have 
studied their effect on the modulus and on the phase of 
the transfer functions. We have found that the uncer- 
tainty on the capacitive actuator gain (6gca) ^^e only 
one that has a relevant impact, whilst the uncertainties 
on the other parameters have a negligible effect. 



A. The gain of the capacitive actuator 

In this section we analyze the effect of the uncertainty 
of the capacitive actuator gain (&Gca)- To this end, for 
a specific value of bccA^ compute the absolute error 
of the modulus (Ho*^) of the system transfer functions 

across the measurement 



and the phase differences (H^ ) 
bandwidth (ImHz to 30mHz): 



V. AN ANALYSIS OF THE UNCERTAINTIES 

The Experiment Master Plan of the mission is aimed 
at determining the physical parameters of the instru- 
ment, characterizing in this way the matrix elements of 
H. These transfer functions depend on several physical 
parameters. Amongst them we mention the stiffnesses of 
the test masses (wi and 0J2, where the subindexes 1 and 
2 refer, respectively, to test mass 1 and 2), the actuator 
gains (namely, the gain of the FEEP actuator, Gfeep, 
and that of the capacitive actuator, Gca), and the in- 
terferometer cross-coupling (i5i2)- In the end, this results 
in a complete characterization of the main four blocks of 
the instrument. Detailed information on the Experiment 
Master Plan and on the accuracy of the estimates can be 
found in Refs. [idmil]. 

Nevertheless, for the calculations presented here it is 
important to realize that some of the parameters of the 
model may be poorly determined or have sizable uncer- 
tainties. Therefore, in our analysis we introduce uncer- 
tainties in each of the most relevant parameters of the 
mission. These uncertainties are represented as 5, and 
the subscript "NOM" stands for the nominal value of 
the parameter: 



He*^ = |H'^|-|H| (23) 
H^'^ = ZH'^-ZH (24) 

where Z stands for the matrix operator that calculates 
the phase of each of the elements of the matrix. In these 
expressions, the superscript "b" indicates that the spe- 
cific transfer function has been calculated with a non- 
zero value of the parameter uncertainty. On the contrary, 
functions without superscript have been calculated with 
the nominal values of all the system parameters. There- 
fore, HJ? calculates the absolute error of the modulus of 
each of the nine functions with respect to its nominal 
value for one specific value of uncertainty (6), and H^,*^ 
gives account for their phase differences. These two ma- 
trices give a quantitative assessment of the error of the 
model due to the uncertainties across the entire measure- 
ment bandwidth. Moreover, the most relevant contribu- 
tion in the error of the model will be due to the error in 
the diagonal terms of the matrix. Therefore, we analyze 
mainly the effects on Hp^^xi2i Hn^^jji and H^^^^-^. 

Fig.|3]displays the results of this sensitivity analysis for 
Hp^^xi2j i-G the first element of matrices Ho'' and H^'' 
for several values of the uncertainty in the capacitive ac- 
tuator gains, ranging from —0.2 to 0.2. The behavior as a 
function of the frequency of the other two elements of the 
diagonal are very similar. In the top panel of this figure 
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FIG. 3: Top panel: error of the modulus of the transfer func- 



tion H% 



with respect to its nominal behavior. This 



frequency-dependent relative error is plotted for different ca- 
pacitive actuator gain uncertainties ranging from —0.2 to 0.2. 
Bottom panel: phase differences in the Hp^_^^^^ transfer func- 
tion for different uncertainties of the gain of the capacitive ac- 
tuator. The phase differences are also calculated for different 
relative gain uncertainties ranging from —0.2 to 0.2. 



it can be seen that the error of the modulus is especially 
relevant below 1 mHz, where the differences in amplitude 
increase up to 48% for 0.6 mHz, when the capacitive gain 
is 0.8 (instead of 1). The changes in modulus are also rele- 
vant between 1 mHz and 7 mHz. In the bottom panel, we 
examine the differences in the phase of the same transfer 
functions. It can be seen that there exist phase shifts of 
15° for a capacitive actuator gain of 1.2 at a frequency of 
1 mHz. These phase shifts are relevant between 0.4 mHz 
and 4 mHz. Such differences produce important biases 
in the estimates of the magnetic parameters. Moreover, 
the effect depends on the excitation frequency. Thus, 
the choice of the right excitation frequency (wq) is a cru- 
cial aspect in the experiment design. We postpone this 
analysis to Sect. |VI[ where we will study which is the op- 
timal excitation frequency. Finally, we also mention that 
similar analyses for the rest of the uncertainties on the 
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FIG. 4: Signal-to-noise ratio as a function of frequency for 
each of the relevant signals of the magnetic experiment, the 
differential channel, X12 — solid line — the rotation about the 
y-axis, 771 — dashed-dotted line — and the rotation about the 
z-axis, 01 — dotted line. 



nominal parameters have been performed, but are not 
shown here for the sake of conciseness. 



VI. THE OPTIMAL FREQUENCY 

Finding the optimal frequency of the sinusoidal cur- 
rents injected in the coils to obtain the magnetic param- 
eters is a crucial issue of the experiment. Actually, as 
it will be shown below, the optimal frequency can be 
obtained from a trade-off between the frequency range 
where the instrument presents a maximum of the signal- 
to-noise ratio (SNR) and the frequency range where the 
instrument is less sensitive to the uncertainties of the ca- 
pacitive actuator gain — see Sect |V A[ 

The SNR across the instrument measurement band- 
width for each of the channels — 0x^2 1 , and o^-^ — is 
shown in Fig. [4] The SNR reaches its maximum between 
0.5 to 1.5 mHz for the displacement reading, and from 
1 to 2niHz for the rotation channels. This is the most 
sensitive band of the instrument. This is confirmed by 
inspecting Fig. [5] where we show the response of the sys- 
tem to the excitation by 4 different sinusoidal currents. 
All these sinusoidal currents have the same amplitude, 
1mA, but they oscillate respectively at 0.5 mHz, 2 mHz, 
5 mHz and 9 mHz. In the top panel, we show the read- 
ings of the differential displacement channel to this set 
of four sinusoids. When exciting at 0.5 mHz the ampli- 
tude is ^ 40 nm, whereas at 2 mHz drops to ^ 5 nm. 
Finally, when the frequency is 5 mHz the amplitude of 
the excursion is only ^ 1 nm. This same effect is ob- 
served in the bottom panel, where we show the Fourier 
analysis of these time series. As can be seen, each of 
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FIG. 5: Top panel: time series of the responses detected at the 
differential channel of the interferometer when we inject 4 dif- 
ferent sinusoidal signals in the onboard coils. The amplitude 
of these sinusoids in all the cases are 1 mA and the frequencies 
are respectively 0.5 mHz, 2 mHz, 5 mHz and 9 mHz. Bottom 
panel: Fourier analysis of the time series displayed in the top 
panel of this figure. 



the readings has a frequency component at loq and a 
second one at 2ujo, as expected. Note as well that the 
2wo components are highly attenuated with respect to 
the main component because they are located at higher 
frequencies. This simple analysis seems to indicate that 
the excitation frequency should be chosen around 1 mHz. 
However, this range of frequencies is where the uncer- 
tainty of the capacitive actuator has the largest impact 
on the estimates of the magnetic parameters — see Fig. [3] 



and section V A Thus, the determination of the optimal 
excitation frequency should be the result of a joint op- 
timization procedure, taking into account both the fre- 
quency dependence of the SNR and the uncertainties in 
the gain of the capacitive actuator. 

To find the optimum excitation frequency we compute 
the estimation error of each of the magnetic parameters 
for different uncertainties of the gain of the capacitive 
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FIG. 6: Global error functions of each of the magnetic param- 
eters with respect to the excitation frequency. The dashed 
black line is the global error function for Mj, — that is, 
EM^ifo) — the black solid line is EMyifo), and the dotted 
black line corresponds to Em^ (fo)- Finally, the solid gray line 
corresponds to E^^{fo) and the dashed gray line to E^^{fo). 



actuator ranging from —0.2 to 0.2. We do this for differ- 
ent excitation frequencies across the entire measurement 
bandwidth. Thus, for each magnetic parameter, we com- 
pute an error function for each gain uncertainty, eb{u!o)- 
Then we add quadratically each of these functions with 
their appropriate weight factor: 



E{ujq) = (\eb{uJo) 



(25) 



where b is the uncertainty in the capacitive actuator gain, 
in percentage. In this way we compute a global error 
function for each of the magnetic parameters, i?A/j.(wo), 
EMyi^^o), EmA^o), E^^loo), and E^^ujo). The absolute 
minima of these functions correspond to the best excita- 
tion frequencies for each of the parameters. 

The global error functions computed in this way are 
shown in Fig. [6] for frequencies from 0.1 mHz to 12 mHz. 
For the case of remanent magnetic moment the er- 
ror function presents a very broad minimum between 
^ 5 mHz and ~ 11 mHz, being the absolute minimum 
at ^ 10 mHz. Note that at lower frequencies the global 
error function grows very abruptly. This occurs because, 
although the SNR of the experiment is larger at these 
frequencies, they are also very sensitive to the biases in- 
troduced by the actuator uncertainty. Note as well that 
the error functions have local minima at around 1 mHz, 
and also a local maximum between ~ 1.2 and ~ 2 mHz, 
following the sensitivity curve of the capacitive actuator 
— see Fig. [3j On the other hand, the optimal frequency 
needed to estimate Xo and Xo lies between 5 and 7 mHz. 
This is because these last two parameters are estimated 
with the 2ujQ component of the Xi2 output, and higher 
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FIG. 7: Statistical distributions of the estimates for the 5 magnetic parameters of the test masses for 3 different simulations. 
The left column shows the statistical distribution of these parameters when the capacitive actuator gain has no uncertainty. 
This simulation is done for an excitation frequency of 5 mHz. The second column shows the results when the uncertainties 
of the capacitive actuator gains of the principal axes are modeled with a normal distribution of zero mean and 0.01 standard 
deviation. This simulation is done for an excitation frequency of 5 mHz, too. Finally, the third column, shows the parameter 
estimation results for the same experiment as in the second column, but for 1 mHz. The x-axis of each subplot shows the 
relative error in the parameter (in percentage). 



frequencies are penalized by the larger attenuation on 
this frequency component. Finally, it is worth mention- 
ing that the phase shift shown in Fig. |3] around 1 niHz 
penalizes the estimation at low frequencies, because the 
components at 2ujo suffer a different and unknown shift 
with respect to the luq component. In summary, the best 
choice of excitation frequency is 5 mHz to estimate Xo 
and Xc and 10 mHz to estimate the three components of 
the magnetic moment. Nevertheless, if only one inflight 
experiment could be performed due to planning restric- 
tions of the mission, the best frequency would be 5 mHz. 
This value is the result of minimizing the quadratic sum 
of the five error functions of the five parameters. 



VII. ROBUSTNESS OF THE ESTIMATES 

Finally, to conclude with our analysis we have stud- 
ied the robustness of our findings. Specifically, we have 
tested the performance of our estimation algorithm un- 
der several circumstances. In order to model statistically 
its performance, we have estimated the magnetic parame- 
ters for 50 different simulated experiments and calculated 
the statistical distribution of the relative errors of each 
parameter. For example, for the case of the relative 



error is computed as: 

where AI^ is the estimated parameter and represents 
its true value. 

Here we present the results of three different simula- 
tions. In the first simulation we excite the coils with 
a 5 mHz sinusoid and we consider that the gain of the 
capacitive actuator is the nominal one. In the second 
simulation we maintain the 5 mHz excitation frequency, 
but in this case the gains of the capacitive actuators of 
the three main axis are modeled with a normal random 
distribution of zero mean and of 0.01 standard deviation. 
Finally, the third simulation is only performed for illus- 
trative purposes. We maintain a random distribution of 
the uncertainty of the gain of the capacitive actuator, 
but we excite the coils at 1 mHz. Note that the analy- 
sis previously explained in Sect. |IV| concluded that this 
frequency should not be used. Consequently, this case 
clearly illustrates the effect of choosing a wrong excita- 
tion frequency. 

The error distributions for each of the simulations pre- 
viously described and for each of the magnetic parame- 
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TABLE I: Standard deviations of the estimated parameters 
for different estimation scenarios. For each of the different sce- 
narios we calculate the ratio between the actual performance 
and the optimal Cramer Rao lower bound. 



Run 




AM^ 


AMy 


AM, 


Axo 


Axc 




CR bound 


0.019% 


0.046% 


0.139% 


0.083% 


0.263% 


1 


(7 


0.028% 


0.067% 


0.156% 


0.084% 


0.557% 


CR ratio 


1.47 


1.45 


1.12 


1.01 


2.11 


2 


a 


0.123% 


0.132% 


0.162% 


0.176% 


0.632% 




CR ratio 


6.47 


2.86 


1.17 


2.12 


2.40 




(7 


0.331% 


0.215% 


0.445% 


0.553% 


6.557% 


3 


CR ratio 


17.42 


4.67 


3.20 


6.66 


24.93 



x-component of the remanent magnetic moment we may 
write: 



N 

1=1 



(27) 



where N is the total number of frequencies used, uji is 
the corresponding excitation frequency, M^^. is the esti- 
mate of Mx at uji and Mx is the final combined estimate. 
In this equation EM^{^j->i) are the weighting factors of 
Eq. ( 25 ) adequately normalized: 



N 



ters are displayed in Fig. [7j and their respective standard 
deviations are listed in Table |T] For consistency, these re- 
sults are checked against the Cramer-Rao lower bound, 
which gives a lower limit for the variance of the estimated 
parameters [TB] . The Cramer-Rao bounds for each of the 
estimates of the magnetic parameters are listed the first 
row of Table |lj Moreover, for each of the simulations pre- 
sented here we also compute the ratio of the standard 
deviation to the Cramer Rao lower bound. For the first 
simulation we obtain variances close to the Cramer-Rao 
lower bound, as expected due to the large SNR. In the 
second numerical experiment we obtain standard devia- 
tions smaller than 0.18% for all the magnetic parameters, 
except for Xc, which is the one with the lowest SNR. In 
this experiment, we are still close to the optimal Cramer 
Rao bound because we minimize the effect of the capaci- 
tive actuator uncertainty. Finally for the third simulation 
we obtain an important degradation of the performance 
of the parameter estimation procedure. In particular, the 
standard deviations are increased by more than 1 order 
of magnitude. The ratio with respect to the Cramer-Rao 
lower bound is also clearly much larger. Particularly, the 
performance of the estimate of Xc is totally unacceptable 
for this experiment, obtaining an estimation performance 
25 times worse than the optimal one. Finally, comparing 
the second and third columns of Fig. [7] — and the sec- 
ond and third sections of Table [J — we confirm that our 
estimation procedure delivers better results (and close to 
optimal) for an excitation frequency of 5 niHz than for 
1 mHz, which was the frequency adopted in the prelimi- 
nary design of the experiment. This clearly demonstrates 
the importance of choosing the appropriate excitation 
frequency. 

However, this is not the most robust estimate that can 
be obtained. In particular, we suggest to use a multi- 
frequency estimation technique, where the properties of 
the test masses are computed using the results obtained 
at different frequencies. In this way the effects of spuri- 
ous or non-modeled effects at a given specific frequency 
can be minimized. This can be done weighting the re- 
sults obtained for each of the magnetic parameters at 
each frequency by the inverse of the corresponding total 
error function given by Eq. (25 1. For instance, for the 



(28) 



This estimation procedure provides an estimate of the 
magnetic characteristics of the test masses that takes into 
account all the limiting factors of the LTP instrument, 
and also delivers estimations which are robust to other 
unexpected (and not modeled) frequency dependent ef- 
fects. 



VIII. MAGNETIC CONTRIBUTION TO 
PROOF-MASS ACCELERATION NOISE 

In previous sections we have shown that the magnetic 
properties of the test masses can be characterized in flight 
with accuracies below 1% for each of the parameters. The 
ultimate goal of measuring these properties is the deter- 
mination of the magnetic contribution to the proof-mass 
acceleration noise. To do this, the estimated values of 
the magnetic field and magnetic field gradient at the po- 
sitions of the test masses are needed. It has been con- 
cluded [T71[TH] that, in the worst case, the magnetic field 
and its gradient can be estimated with an accuracy bet- 
ter than 10%. The spatial nonhomogeneities of the field 
are included in this uncertainty. Thus, an important goal 
consists in assessing how these uncertainties project into 
the precision of our estimate of the magnetic contribution 
to the total differential acceleration reading in LPF. Ac- 
cordingly, we propagate the errors of the magnetic char- 
acteristics and the magnetic field and gradient into the 
calculation of the magnetic acceleration noise: 



Ctotal 



a/ 



TV 



^ V ds,. 



(29) 



where / is the total magnetic contribution to the ac- 
celeration noise and Si are the several sources of error. 
In our case, the total magnetic acceleration noise is the 
sum of the fluctuation of the magnetic field and magnetic 
field gradient of the spacecraft, the down converted AC 
magnetic fields, the interplanetary magnetic field fluctu- 
ations, and the Lorentz force contributions jl9j . whereas 



11 



TABLE II: Budget of the contribution of the magnetic field effects to the total acceleration noise and their uncertainties. 
These uncertainties are computed using the expected error of the magnetic characteristics reported in this paper and using the 
expected error of the magnetic field determination reported in Refs. [171 118j 



Contribution 


Differential acceleration noise [m s ^ Hz ^^^] 


Fluctuation of the spacecraft magnetic field 


(0.680 ±0.096) xlQ-''^ 


Fluctuation of the spacecraft magnetic field gradient 


(1.097 ±0.108) xlO"^^ 


Down converted AC magnetic fields 


(1.265 ±0.254) xl0~^^ 


Interplanetary magnetic field fiuctuation 


(1.701 ±0.241) xlO"^^ 


Lorentz force 


(0.013 ±0.001) xlO"^^ 


Total 


(2.775 ±0.425) xlO"''' 


Requirement 


12.0 xlO"'^ 



the sources of error are the uncertainties of the magnetic 
field and its gradient and those of the remanent magnetic 
moment and susceptibility of the test masses. The results 



obtained using Eq. ( 29 ) and the computed uncertainties 
are shown in Table |II[ The total magnetic contribution, 
which is 2.775 x 10"^ m s^^ w-^-i/s ,c ^o+oi-mino^ 



Hz is determined with 

15 ^ £,-2 



with an accuracy of 0.425 x 10~^^ m s~'^ Hz~^/^. This 
means that the magnetic contribution to the total noise 
can be estimated with a fair accuracy, and therefore, can 
be subtracted from the main acceleration reading, with 
a relative error 15% across the whole LTP measurement 
bandwidth. This represents an enhancement of one order 
of magnitude with respect to previous studies. 



IX. SUMMARY AND CONCLUSIONS 

In this paper we have studied how the magnetic charac- 
teristics of the test masses onboard LISA Pathfinder can 
be determined. This is essential to estimate the magnetic 
noise contribution to the entire noise budget and, most 
importantly, to subtract this noise from the displacement 
reading. The estimation of M^, My, M^, Xo and Xe is 
done by injecting a controlled magnetic field at the po- 
sition of the test masses. The field is generated by a 
sinusoidal current circulating through the two onboard 
induction coils placed at each side of both test masses. 
The induced magnetic field results in magnetic forces and 
torques on the test masses that excite their dynamics. 
We have shown that the force acting on the test masses 
has two frequencies, while the torques oscillate at single 
frequency, allowing to estimate the properties of the test 
masses by an adequate processing of three of the readings 
delivered by the interferometer. These readings are the 
differential displacement of both test masses {0x12)^ the 
rotation of test mass 1 about the y-axis (o,,^) and that 
about the z-axis (o^^ ) . We have also shown that the time 
series received from the satellite's telemetry need to be 
whitened and split to minimize the low-frequency effects 
inherent in the operation of the instrument. This way. 



the magnetic parameters can be estimated by a classi- 
cal single-channel least-square technique once the effects 
produced by cross-talks are determined and corrected. 
Additionally, we have assessed the sensitivity of the es- 
timation procedure to the uncertainty in the gain of the 
capacitive actuator of the instrument. This effect showed 
to be very relevant and, most importantly, it has been 
found to depend on the excitation frequency. Moreover, 
the SNR of the received signals also depends on the fre- 
quency of the injected signal. Accordingly, we have also 
presented a joint optimization analysis that takes into 
account these two factors, leading to the conclusion that 
the optimal excitation frequency for a joint experiment 
is 5mHz. Performing the experiment at this frequency 
allows to estimate the magnetic characteristics without 
being affected by the likely uncertainty in the capacitive 
actuator gain. In this case we obtain parameter vari- 
ances smaller than ~ 0.7% when the deviations of the 
gain of the capacitive actuator are ~ 1%. Using all the 
previously explained steps and adopting this excitation 
frequency, the estimation turns out to be more accurate 
than that obtained using the preliminary design of the 
experiment, for which a frequency of 1 mHz was adopted. 
Moreover, we have suggested that a multi-frequency es- 
timation technique could deliver estimates of the highest 
quality, enhancing the robustness of the experiment in 
front of non-modeled frequency-dependent effects. Fi- 
nally, taking into account the aforementioned results, to- 
gether with the results about magnetic field estimation 
presented elsewhere, we estimate that the total magnetic 
contribution to the proof-mass acceleration noise can be 
determined with a ^ 15% error level across the whole 
LTP frequency band. 
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